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Tunneling of two particles in synchronous and asynchronous regimes is studied in the framework 
of dissipative quantum tunneling. The critical temperature T c corresponding to a bifurcation of 
the underbarrier trajectory is determined. The effect of a heat bath local mode on the probability 
of two-dimensional tunneling transfer is also investigated. At certain values of the parameters, 
the degeneracy of antiparallel tunneling trajectories is important. Thus, four, six, twelve, etc., 
pairs of the trajectories should be taken into account (a cascade of bifurcations). For the parallel 
particle tunneling the bifurcation resembles phase transition of a first kind, while for the antiparallel 
transfer it behaves as second order phase transition. The proposed theory allows for the explanation 
of experimental data on quantum fluctuations in two-proton tunneling in porphyrins near the critical 
temperature. 
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I. INTRODUCTION 

Quantum tunneling dynamics of a particle in- 
teracting with a heat bath is one of the im- 
portant problems of modern condensed matter 
physics 1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 



The interest to this problem is related to studies in 
low-temperature superconductive tunnel junctions^ 4 ^ 
a dissipative quantum tunneling in crystals^ and low- 
temperature chemical reactions iS 4 *2Si2Si2L22i22i 3 A In low 
dimensional systems an effective mass approximation 
often fails. Thus, quantum tunneling with dissipation 
becomes an important tool in the description of elec- 
tron transferrin In many physical processes tunneling 
of two particles occurs. For example, Semenov and 
DakhnovskiiSSiS 9 . found the bistability of the tunneling 
trajectories in a two-proton transfer. Later, Benderskii 
and co-workers extended their investigations to different 
two-dimensional potentials. 30 Some features of the 
two-dimensional tunneling dynamics were studied in 
interacting Josephson junctions,^ 

Another example of two-particle tunneling is a two- 
proton transfer in porphyrin systems ^LSSiS 9 . Porphyrins 
are important molecules in Biology^ and a new area of 
Electronics - molecular wires and devices^ There are 
experimental data^S^L^ 3 . that clearly indicate that a 



tunneling instability occurs at some critical temperature. 
Such a behavior in the rate constant demonstrates the 
existence of bifurcations in two-dimensional underbarrier 
trajectories. Apparently, this effect requires a thorough 
investigation when the protons interact with a thermal 

2f> 29.30 



In this paper, we continue to study bifurcation ef- 
fects in two dimensional tunneling. In particular, we 
consider two-proton correlations of different types within 
the framework of dissipative quantum tunneling (instan- 
ton) approach. Moreover, we discuss a fine structure of 
bifurcations for systems with various potential energy 
surfacesi 27 i 28 i 29 i 30 i 33 i 34 In Section II, we introduce two- 
dimensional model potential energy surfaces for a pair of 
interacting particles. In Section III, we study the effect 
of temperature on the tunneling rate. In Sections IV and 
V, we calculate the rate for the parallel and antiparallel 
tunneling and provide the analysis of the origin of bifur- 
cation. The effect of a promoting mode (an environment) 
is studied in Section VI. 
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II. TWO-DIMENSIONAL POTENTIAL 
ENERGY SURFACES 



Consider two charges that tunnel in two indepen- 
dent double-well potentials Ufa) and Ufa) presented 
as follows. 7 i 27 i 28 i 29 



Ufa)= l -uu 2 fa+a) 2 9{-q l 



~AI+-Lu 2 fa-b) 2 



Ofa 



i = l,2, 



(1) 



where the sum a + b determines the length of a " link" in 
the corresponding macrocluster fragment; AI — Lu 2 (b 2 — 
a 2 )/2 is a bias (an asymmetry parameter of the poten- 
tial); Ofa) is a step function, and lu is the frequency (see 
discussion in Ref. l27l) . The mass is absorbed into the 
definition of q. 

The interaction between two charges, e.g., protons, is 
considered in a dipole-dipole approximation 2 ^ 




Viatfa,q2) = -77(91 - Q2) 2 , 



(2) 



where a is a positive constant. We use the same interac- 
tion potential as in Ref. |y. 

Thus, the total two-dimensional potential energy sur- 
face for parallel tunneling normalized by lu 2 , is given by 



FIG. 1: Asymmetric potential energy surface Q for parallel 
tunneling; a = 2, b = 2.5, a* = 0.0001. A and B indicate 
the initial and final states of the particles, respectively. The 
minimum of the potential at B is lower than that of at A. 
The other two (intermediate) minima are lower than those of 
at A and higher than at B. 



U p fa,q 2 ) 



2U p fa,q 2 ) 



= fa + a) 2 9(- gi ) 



+ [-(b 2 - a 2 ) + fa - b) 2 } 9 fa) + fa + a) 2 9(-q 2 ) 
+ [-(b 2 - a 2 ) + fa - b) 2 } 9 fa) - ^fa ~ q 2 f. 



(3) 



Here, a* — 2a/ lu 2 is the dimensionless parameter, a* < 
1, a ~ e 2 /(ei?o)> e is the electron charge, Rq is the 
separation distance between the reaction coordinates q\ 
and qi of the tunneling particles; e is the dielectric con- 
stant. The form of potential energy surface is shown 
in Fig. 1. 

For antiparallel transfer, the two-dimensional potential 
energy with the interaction term can be defined as 



£4(91,92) 



214(91,92) 



fa+a) 2 9{- qi ) 



<b 2 



2 ) + fa-b) 2 ] 9fa) + fa-a) 2 9fa) 



+ [-(b 2 -a 2 ) + fa + b) 2 ]9(-q 2 ) 



-fa-q 2 ) 2 1 (4) 



where a* = 2a /lu 2 is a dimensionless parameter, a* < 
1. The potential J3J is depicted in Fig. 2. The main 
difference between the potential energy surfaces (3) and 
(4) is in the initial location of the second particle, ±a. 

The potential energy 10 can be referred to as "paral- 
lel" while the potential energy (0} as " antiparallel" . We 
also define a "symmetric" potential energy as a particu- 
lar case of the potential (0} under the condition a = b, 



i.e., 



U s fa,q 2 ) = 



2*4(91,92) 



= fa+a) 2 9(- qi ) 



+ (9i - a) 2 9fa) + fa - a) 2 9fa) + fa + a) 2 9(-q 2 ) 

^-(91 - 92) , 



(•5) 



where the constant a** < 1. The potential energy (J3J) is 
presented in Fig. 3. 

As shown in Refs. 28 29 33, such model potential en- 
ergy surfaces well describe the dynamics in porphyrins. 

In many practical cases, the effect of a bath on parti- 
cles' tunneling should be also included. In the next sec- 
tion we consider two particles embedded into a harmonic 
medium and linearly interacting with the bath modes. 



III. TWO-PARTICLE TRANSITION 
PROBABILITY 

We assume that two particles independently interact 
with a harmonic bath. Such interaction is considered in 
a bilinear approximation. The dynamics of the environ- 
ment is described by the oscillator Hamiltonian (we use 
H = 1, ks = 1 units with the oscillators' masses equal to 
one), 



p 



?)• 



(6) 



3 







FIG. 2: Asymmetric potential energy surface £IJ f° r antipar- 
allel tunneling; (a) a — 2, b — 2.3, a* =0.1 (top panel); (b) 
a = 2, b — 2.3, q* = 0.5 (bottom panel). A and B indicate 
initial and final states of the particles. The minimum of the 
potential at B is lower than that of at A. The other two 
(intermediate) minima are higher than those of at A and B. 



FIG. 3: Symmetric potential energy surface JjjJ; (a) a — 2, 
6 = 2, &** = 0.1 (top panel); (b) a = 2, 6 = 2, a" = 0.5 
(bottom panel). A and B denote the initial and final states 
of the particles. The minimum of the potential energy at B is 
equal to that of at A. The other two (intermediate) minima 
are higher than those of at A and B. 



Each of the tunneling particles (electrons or effective 
charges) interacts with the oscillator bath in the following 
way: 

Qi) = 9i E C *Q*> ^pihfe, Qi) = Q2 ]T CiQi, 

i i 

(7) 

As in Ref. l28l we are interested in the transition proba- 
bility per unit time or, more precisely, only in its expo- 



nential part which can be written in the Langer's form, 

ImZ 



r = 2T- 



ReZ 



(8) 



In such a consideration, metastable levels can be pre- 
sented as 



r=-2Im£, E = E -iT/2. 



(9) 



4 



It should be emphasized that Eq. (8) is valid only for 
temperatures below the crossover temperature while for 
higher temperatures a different prefactor is required^ 3 - 

Equation iJSJ is obtained by gener- 
alizing the expression @ to nonzero 
temperatureaWiWiVi^iVV'VV'VfiiiT 

2J2e- Em/T lmE i _ 2TImJ>-' B */ T _ 2TlmZ 



-Eoi/T 



Re£, 



-Ei/T 



ReZ 



(10) 

Here, i labels the energy levels in the metastable state, 
Z is the partition function of the system, and T is the 
temperature. 

To calculate T, it is convenient to present Z in the form 
of a path integral 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 . 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 

Z = Tlf DqiDq 2 DQ t exp[-5{gi, q 2 , Qi}}. (11) 

i 

Here, S denotes an underbarrier action for the total sys- 
tem. The imaginary part ImZ is due to the decay of the 
energy states in the initial well. The validity of this ap- 
proximation requires dissipation to be strong enough so 
that only an incoherent decay occurs ^ The imaginary 
part in the partition function with a double well potential 
energy can be also explained due to the strong dissipa- 
tion to the bath. Indeed, the particles do not come back 
to their initial state. Coherent oscillations can happen 
only if the interaction with bosons is weak enough 8 or 
the bath is in noncquilibrium state^ 5 ^ 7 - 

The integral (|llf> can be performed over phonon 
coordinates 2 ^, resulting in 



1 -2 

2 91 



1 -2 

2<?2 



V(qi,q 2 ) 



S{qi,q 2 } = J dr 

-13/2 

13/2 

+ J dr'D(r-T'){q 1 (T)+q 2 (T)][q l (r , )+q 2 (T')}) 



-13/2 

where 



D ^ = l E D ^y 



,(12) 



(13) 



n— — oo 



(3 = h/(ksT) is an inverse temperature (below we assume 
that h = 1 and ks — 1), v n = 2irn/ fi is the Matsubara's 
frequency, and 



D{v n ) 



E: 



(14) 



A trajectory that minimizes the Euclidean action S can 
be found from the equations of motion. In particular, we 
embark on the antiparallel tunneling 

0/2 

- q\ + Q 2 oqi + a x q 2 + / dr'K{T - r')[gi(r') + q 2 {r')] 



-P/2 

+LU 2 a9(- qi ) - u 2 b6( qi ) = 0, (15) 



0/2 

- q 2 + fio?2 + &iq 2 + J dr' K (t - T')[gi(r') + q 2 {r') 



-0/2 



-LU 2 a8(q 2 ) + LU 2 b9(-q 2 ) = 0. (16) 
In Eqs. I|15l) and l|l(jfl the kernel K is defined by 



-. oo 

= p E ^ 



(17) 



Here, £ n is the Green's function defined by Eq. (14) with- 
out the zero-frequency term, 



£>("») = -E-?+& 



(18) 



Thus, we seek for solutions of Eqs. (|T5|) and l |T5jl by 

expanding the trajectories qi{t) and q 2 (t) into Fourier 
series, 



1 oo oo 

ft = \ E tf } ^ T . ^ = \ E 



/3, 



^e^. (19) 



Introducing the renormalized frequency and interaction 
constant, 

respectively, and substituting Eqs. ljT§|l into Eqs. i(T5|) 
and Ijl6|l . we obtain that for n = 0, 

(l) . (2) 2w 2 (a + &)e 

Si + OLi 

(1) (2) _ 2^ 2 a/3 4^ 2 (a + 6)r 

% % — 752 ~ 1 f^2 - ' v zl ; 



and for ti^O, 



2) _ 2a; 2 (a + &)(sinz/ n Ti — sin ^„r 2 ) 



>(^ 2 + ft 2 + 6^ + 2£„) ' 



(l) _ (2) _ 2a; 2 (a + 6) (sin tyri + sin ^»t 2 ) 

Here, we have introduced the following notation: 

e = n-T 2 , t = (n +t 2 )/2. (23) 

The time instants Ti and r 2 , at which the particles pass 
the top points of the barrier, are determined from the 
following equations: 



?i(n) = 0, <? 2 (t 2 ) = 0. 



(24) 



Equations l|24|) allow us to change the argument of the 9- 
function. Namely, instead of the q\- and ^-dependencies, 
we can use a time-dependent ^-function. This reduces 
Eqs. (|15|l and (|16|) to a linear form. 
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Finally, substituting the trajectory determined from 
Eqs. O, J2U and (E3> mto Eq. JT^J), we arrive at the 
following expression for the instanton action: 



S = 



4w 4 a(a + 6)r uj 4 (a + b) 2 e 2 4w 4 (a + 6) 2 t 2 



O 2 



ai 



/3(ng + a a 



/?(tt 2 - 5 a 



8cj 4 (a + 6) 

is" 



2 °° 

E 

n=l 



sin 2 v n T ■ cos 2 (v n e/2) 



("2 + - 
cos 2 ^„r • sin 2 (t' n e/2) 



(^ +0 2 + a 1 + 2£„K 2 



.(25) 



IV. 



PARALLEL PARTICLE TUNNELING 



For the case of parallel particle tunneling, the Eu- 
clidean action S can be determined similarly to antipar- 
allel tunneling [see Eqs. i|15fl and (|16l) ]. The trajectory 
minimizing the Euclidean action (instanton), can be de- 
termined from the equations of motion. As in the pre- 
vious section, we seek for solutions of these equations in 
the form of the Fourier expansion Ijl9(l . The time instants 
n and r 2 are determined by Eqs. J23J- 

In the case of parallel tunneling particles [the potential 
energy ©], the resulting Euclidean action is given as 
follows: 

S = 2a(a + b)(n + t 2 )lu 2 - -L 2 (a + 6) 2 (ti + t 2 ) 2 

P 



tu\a + 6) 2 (n - r 2 ) 2 



2a)(3 



2uj A {a + bf 




E 

n=l 



(sin v n T\ + sin v n T2) 2 
(sini^ri - smv n T 2 ) 21 



(26) 



where £ n is defined by Eq. I|18|l . 

Below, we use the following notation: 

£ = e*u! = (n — T2)u), 



t = 2t*uj = (n + t 2 )uj, 
P* = Puj/2, 
a* = 2a/uj 2 , 



and assume that 



b* = b/a, 



b> a. 



In the absence of interaction with an oscillator bath, 
i.e., at £„ = 0, the action (|26|1 as a function of the pa- 
rameters e and r yields 



(a + b) 2 uj ( 4ar 



a + b a + b 



1 



1 - a* 



+ (t \ £ \) a +cothj3 * _ sinir 1 /?* [cosh(/3*-r)coshe 
1 — a* 

+cosh(/3*-r) - cosh(/3* - |e|)] 
-(1 - a*)~ 3/2 [-coth(/Vl - a*) 
+sinh -1 (/Vl - a*)[cosh((/3*-r)Vl - a*) 

x(cosh(eVl - a*) - 1) + cosh((/3* - |e|)Vl - a*)]] j. (27) 

As soon as the trajectory is found, Eqs. I|24|) can be 
presented in the following form: 

sinhe[coshrcoth/?* — sinhr — coth/3*] + 



l-a 



■ sinh(eVl-a*) [cosh(rv / l-a*)coth(/3Vl-a*) 



3 - 



-sinh(iVl - a*) + coth(/3*Vl - a*)] = 0, 
4 1 

- + coshe[sinhTCoth/3* — coshr 
— 1] + sinhrcoth/3* — coshr 



1 + 6* l-a 
1 



l-a* 



cosh(e\/l — a*) [ sinh(TvT— a*) 



Xcoth^Vl^a*) - cosh(rVl - a*) + 1] 

— [ sinh(rVl - a*) coth(/3* y/l - a*) 

1 — a* 



-cosh(rVl - a*)} = 0. (28) 

Simple analytic solutions of Eqs. 1|28J) are obtained in the 
particular case when 



(n - t 2 V = 0, V/3, a < J 2 12, 



7i = t-) — — = — arcosh 
" 2d 2oj 



l-b* . Blo 

sinn — 

1 + 6* 2 



P 



(29) 



However, a complete analysis requires numerical studies. 

At sufficiently low temperatures, ujj3 3> 1, for 1 < 
b/a < 3, and 

b-a 2a „ _ 2(6 -a) 
2(6 + a) - w 2 K " c = 36- a ' 

we finally obtain, with the exponential accuracy, 



-Ti/1 — a* 



1 



+ 6* 



1 — a 



xn + (l-a*) 



- n-*^ 1 /(i-v / T :r ^ T )r_. 



1 - a* 



(3- 



1 



1 + 6* 
-l 



1 - a 
1 

l-a* 



-)/(l - y/T=&)] 
1 



1 — a 



- (30) 
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The solution (|30|1 is valid at 

rVT 



f3>(3 c 



a* 



CO 



(31) 



We point out that an approximate solution can be 
found for large values of the parameter b* (and small 
a*). However, below we focus on the more important 
solution Q30|l . The analysis indicates that there are no 
perturbative solutions of Eqs. I|28|) at low temperatures 
and small s. 

For e = [see Eq. fZ5)l]. the action (J37J) results in 



uj(b 2 — a 2 )arcosh 



b* - 1 , w/3' 

sinn— 

b* + 1 2 



V(& 2 



a 2 )/3 + w(& + a) 2 

1/2 



cosh 



w/3 



1 



(b* + l) 2 2 



sinh X T- 



(32) 



The action H32|) coincides (up to a factor of two) with that 
of calculated in Ref. |22l Thus, we have calculated the 
two-particle Euclidean action for the case of synchronous 
parallel motion of the two interacting particles. 

In the symmetric case (b* — 1), the action reduces 

to 



S 



o 

a — b 



. 9 , w/3 
= 4w tanh— . 



(33) 



At b* > 1, the character of the temperature dependence 
is almost the same. 

At e 7^ 0, the corresponding action S can be obtained 
by substituting Eq. JHDJ into Eq. (for brevity, this 

cumbersome expression is omitted). A simple analysis 
shows that S E ^o < S e =Q. Moreover, it appears that the 
difference AS = S £7 io — S £ —o has a maximum at to (3 1. 

The tunneling paths flt^l are presented by the solutions 
P§)l and (fTTTTf> - At the critical point, (3 = (3 C , defined 
by Eq. H31|) . a relatively abrupt change in the dynamics 
results in a splitting of the single trajectory (e = 0) into 
two (e 7^ 0), as shown in Fig. 4. 

At (3 > (3 C , i.e., for temperatures below the critical 
temperature, T < T c , only a split trajectory (e ^ 0) 
occurs since S e ^o < S £= o (see Fig. 5). 

When (3 < f) c , i.e., at T > T c , and a > a c [see 
Eq. iJSDJ], the solution is determined by a single tra- 
jectory. We have also found that the single trajectory 
solution holds in the whole temperature range for a sym- 
metric potential energy (b* = 1). 

The above consideration is based on the analytic so- 
lutions H29|) for Tx,s near the critical point. A complete 
analysis of Eq. Q28f) requires more detail. 

The numerical study of the set of transcendental equa- 
tions 128|) reveals remarkable features of two-dimensional 
tunneling. The value of the critical parameter a c de- 
creases with temperature as shown in Figs. 6, 7, and 8. 
The temperature dependence of a c indicates to the exis- 
tence of a finite low-temperature limit. It becomes clear 




FIG. 4: Trajectories (a single path is characterized by s = 
and a splitted path is characterized by e 7^ 0) for two parallel 
moving particles, at uof3 ^> 1; (l)-(4) label the projections of 
the minima of potential energy 



2.5 



1 . 5 




0.2 



FIG. 5: Instanton action as a function of the interaction 
parameter a* = 2a/oj 2 of the two parallel moving particles, 
at uj(3 ^> 1; S p = S/{ua 2 ); (1) is the single trajectory; (2) is 
the split trajectory. 



that the weaker the interaction between the tunneling 
particle, the lower temperatures are required for the syn- 
chronous tunneling. Since the synchronous tunneling is 
valid only for T < T c [see Eq. (|31H ] the dependence of (3 C 
on the interaction parameter reveals that the minimum 
occurs at a* = 0.2. The nonmonotonic behavior of (3 C as 
a function of a* can be explained by the a* -dependence 
of t as well. 

Additionally, the curve in Fig. 8 exhibits an anomalous 
behavior (an increasing part of the curve) that can be ex- 
plained in the following way: a two-dimensional potential 
energy for parallel tunneling is evidently deformed by an 
increase of the interaction constant. Indeed, the minima 
of the potential surface become lower and the distance 
between them is larger. Small deformations in the po- 
tential energy surface change the contribution produced 
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FIG. 8: Parameter /3 C is presented as a function of the inter- 
action parameter of the two tunneling particles. 
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FIG. 6: Numerical solutions of transcendental equations 1281 . 
In addition to the studied solution n = T2, with larger /3 there 
are additional solutions, n 7^ T2, shown in panels (b) and (c). 
Panel (b) reveals a single additional solution corresponding 
to the lower value of the Euclidean action. Panel (c) shows 
two different additional solutions with the one corresponding 
to the lower value of the action. 




FIG. 7: Parameter a c as a function of inverse temperature, at 
the fixed value of the frequency uj and asymmetry parameter 
b*. 



by the temperature increase towards a synchronous char- 
acter of the tunneling. However, for sufficiently large 
deformations of the potential, the situation is a quite 
opposite. The significant deformations are in favor of 
synchronous transfer. Thus, an increase in the interac- 
tion constant provides a similar effect when temperature 
increases. Such a behavior takes place up to the certain 
value of the interaction constant, beyond which the po- 
tential energy surface becomes strongly perturbed. Thus, 
Figs. 7 and 8 are complimentary to each other. Conse- 
quently, Fig. 8 can be viewed as a bifurcation diagram. 
Indeed, the region below the curve corresponds to the 
synchronous tunneling, while the region above the curve 
corresponds to the asynchronous one. 

Fig. 6 demonstrates how the critical parameter (3 C 
changes from synchronous to asynchronous values with 
temperature. Figs. 6b and 6c reveal the existence of the 
additional points, referred to as the bifurcation points, 
corresponding to the change of tunneling regimes. For 
the two bifurcation points shown in Fig. 6c, only one is 
physical because it corresponds to the lower value of the 
action (the second point is metastable). However, for suf- 
ficiently small differences between the values of the action 
at these bifurcation points, both close points contribute, 
thus leading to corresponding "fluctuations" in the sys- 
tem during the change of the regimes. For lower temper- 
atures, such fluctuations become negligible because two 
bifurcation points contributes differently. Thus, a stable 
character of the asynchronous tunnel transfer is achieved 
due to a much lower value of the action at one of the 
bifurcation points. 

For chemical reactions mentioned in Introduction, the 
effect of the change in the tunneling regimes reveals as 
a cleavage in the experimental temperature dependence 
of the tunneling constant. Fine details of this instability 
have not been studied yet. Nevertheless, the very exis- 
tence of instability at the edge of the bifurcation can be 
explained as the result of specific fluctuations. A numer- 
ical analysis of the bifurcations in antiparallel tunneling 
is given in the next section. 
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V. TWO-DIMENSIONAL ANTIPARALLEL 
TUNNELING 

For antiparallel tunneling of two particles [see the po- 
tential energy , the instanton action as a function of 
the parameters e and r is determined by Eq. (|25|l . For 
=0, we obtain 



S = 



LOT{b 2 -a 2 ) uj{a + bf 



1-3* 



e 1- 



1 



1-5 



sinh(|£|\/l — a*) . , , cosh(eVl — a*) + 1 
' (1 o*)3/ 2 - Slnh " £ " + (l-a*)3/2 

x[sinh(/3Vl -a*)] _1 [cosh((/3* - r)\/l - a*) 

-cosh(/?Vl - a*)] 

coshe — 1 



sinh/3* 



[cosh(/3* - r) + cosh/3*] ^ .(34) 



The parameters e and r are found from the following set 
of equations [see Eq. H3J]: 

— sinhe[coth/?* + coshrcoth/3* — sinhr] 



1-3* 



-sinh(eVl - 3*) [coth(/?Vl - 3*) 



sh(rVl - <5*)coth(/3Vl - 3*) 

+sinh(r\/l - 5*)] = 0, 
4 1 



-1 - 



(1 + 6*)(1-S*) 1-3* 
+ (coshe — l)(sinhrcoth/3* — coshr) 

+coshe H — { [cosh(e\/l — 3*) + 11 

x [smhfVVl — 3*)coth(/3*Vl — 3*) 
h(r\/l - 3*)] -cosh(eVl - 3*)| = 0. (35) 



-cos 



Simple analytic solutions of Eqs. (35) can be obtained in 
the following form: 

£ = (n - t 2 )lj = 0, V/3, 3 < u 2 /2, 

T 

T\ = T 2 



arcosh 



- — — sinhl— — v 1 — a*) 
1 + 6* v 2 ; 



2u) 

f.(36) 



Similarly to the case of parallel tunneling, we obtain 
that at low temperatures, lo(3 S> 1, with the exponential 
accuracy, 



A{1 - 3 a 



a/7 



1 - (1 - a*) 1 /'y(A/'y - (1 - 3*)- 1 ) ' 



Ae 



Ty/l — a* 



1-3* 



(37) 



Here. 

A=-l 



(1 + &*)(!- 3*) 1-3* 



1-vT 



a*, 



while 3*, 6*, e, and r are determined as for the parallel 
transfer. 

The solution 1)37(1 is valid at a* cl < 3* < 3* 2 , where 
the lower and upper bounds a* cl and 3* 2 are derived from 
a cumbersome transcendental equation (for brevity it is 
not presented here). Particularly, in the symmetric case, 
b* = 1 , we obtain the condition in a simple analytic form, 
1/4 < 23/cj 2 < 1. 

Furthermore, an approximate solution can be found for 
large values of the parameter b* = b/a (and small 3*). 
However, we restrict our analysis to the more important 
physical solution {57jl. 

The /3-dependent solution l|37|) is valid for (3 > (3 C , 
where 



1 



In- 



Ail-ti*) 1 ^ 



(1-3*)V7(4 



r) 



(38) 



At uj/3 ^> 1, the solutions of Eqs. can be found 
perturbatively (for small e) with given values of the pa- 
rameters (b — a)/(b + a) and 3*. At e — [the solution 
(HSU] , the action J3U yields: 



5-^— ^ arcosh 
(1 — a*) 3 / 2 



b + a 



a . u}(3y/l — 3* 
sinn 



co 2 f3(b 2 - a 2 ) u(b + a) 2 



2(1 -3 2 ) (1-3*) 3 / 2 



coth- 



^Vi 



sinh 



2 

,2\!/2' 



(39) 



For the symmetric potential (a = 6) and s = 0, we 
obtain that (see Fig. 7) 



S 



4uja 2 



(1-3*) 3 / 2 



tanh- 



(40) 



We do not present here a cumbersome expression for 
S £i to which one can obtain by substituting the solutions 
r and e into Eq. (|34|l . A simple analysis reveals that 
S £ ^o > Se=o for P > p c and for relevant 3*. Similarly 
to parallel transfer, the tunnel paths can be found from 
Eqs. l|3l)|) and P7jl. These trajectories on the (91,92)- 
plane are shown in Fig. 9. 

As for parallel tunneling and at (3 > /3 C , the pair- 
tunneling changes from a single to double trajectory 
regime. In contrast to the parallel tunneling, such a 
splitting occurs at any values of the parameters of the 
potential. At [3 > (3 C , we have S e ^o > S £= o and, there- 
fore, S £= q determines the tunneling rate. At (3 < [3 C , the 
two degenerated trajectories are transformed into a sin- 
gle trajectory, qi = —q 2 , corresponding to synchronous 
antiparallel transfer. 

For single particle tunneling, there is only a single tun- 
neling path (instanton) minimizing the action. Hence, 
there are two different types of trajectories for the pair 
of interacting particles. Namely, the main contribution 
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+ <l2 




FIG. 9: Trajectories (the basic path characterized by e = 
and the split one is characterized by e 7^ 0) at u}j3 ^> 1 of two 
antiparallel tunneling particles; (l)-(4) denote the projections 
of the minima of the potential energy surface Ua,(qi,q2) de- 
fined by Eq. I0J. 



to the instanton action is determined either by the sin- 
gle or by the double-degenerated paths depending on the 
value of (3. We also point out that in the case of the par- 
allel tunneling, the particles do not simultaneously pass 
the top points of the barrier, t\ =/= t 2 . This means that 
the tunneling transfer is asynchronous. 

At small values of the interaction parameter a* [see 
Eq. and at temperatures such that (5 < (3 C [see 

Eq. lj2U], there is no splitting of a single path (qi = q 2 ). 
Therefore, the particles pass the top of the barriers at the 
same instants (ri = t 2 ). Consequently, the transfer of the 
particles is synchronous. The temperature dependence 
for the antiparallel transfer action is plotted in Fig. 10 at 
various a*. 

The type of the interaction given by Eqs. ©-©I is 
such that it does not affect the motion along the "cen- 
ter of mass" coordinate, qi — q 2 . For this reason, the 
Euclidean action is independent of the interaction pa- 
rameter as for parallel transfer. Since the state of the 
interacting system characterized by a maximal value of 
the relative coordinate, q\ = —q 2 , is preferable (as it 
provides the lower action), it becomes clear that the in- 
stanton action decreases with the interaction parameter 
in the parallel transfer along the degenerated tunnel tra- 
jectories, and increases with the interaction parameter 
for the antiparallel tunneling. 

For antiparallel tunneling, synchronous transfer (n = 
t 2 ) takes place, while asynchronous transfer is forbidden 
due to the greater contribution to the Euclidean action 
(see Fig. 10). 

The validity condition for weakly interacting 
instanton-antiinstanton pairs, in the adiabatic ap- 
proximation, was discussed in Ref. l27l 




1 2 3 4 5 6 



FIG. 10: Instanton action for antiparallel tunneling (e = 0, 
a = b) as a function of an inverse temperature at the two 
different values of the interaction parameter; 5* a = 5/(4a;a 2 ), 
w = uj/3/4; (1) a* = 0.5; (2) a* = 0.1; (3) a dashed line 
corresponds to the action 1331 1 for the parallel transition (e = 
0, a = b). 

The numerical analysis of transcendental equation l|35(l 
reveals interesting features for a transition region be- 
tween the tunneling regimes, i.e., a fine structure near 
the first bifurcation point for the antiparallel transfer. 
The numerical results are presented in Fig. 11. We found 
that, in addition to the first bifurcation point character- 
ized by the two solutions (Fig. 11a), there exist additional 
bifurcation points at lower temperatures, i.e., the four 
pairs (Fig. lib), the six pairs (Fig. 11c), and there are 
even twelve pairs of additional solutions at (5* — 19.2009 
(a* = 0.05), etc. We refer to this phenomenon as multi- 
plication of bifurcations or a cascade of bifurcations. Such 
an effect resembles a scenario of transition to chaos. 

Although the synchronous regime is preferred due to 
the minimal instanton action, in a certain temperature 
range its value is comparable with those of correspond- 
ing to the cascade solutions. As a result, quantum fluc- 
tuations of a nonregular character occur in contrast to 
the parallel transfer. The antiparallel tunneling is, thus, 
characterized by the instability of the transition due to 
the synchronous to asynchronous behavior. Such insta- 
bilities are similar to continuous second-order phase tran- 
sition, while the parallel tunneling is viewed as a step 
process similar to phase transition of a first order (see 
Fig. 6). The dependencies shown in Figs. 7 and 8, j3 c {a) 
and a c (/3), for the antiparallel transfer are found to be of 
the same character as those of the parallel transfer. 

In summary, our investigation reveals a quite compli- 
cated fine structure of the transition for parallel and an- 
tiparallel tunneling of two particles with different degen- 
erate trajectories leading to the bifurcation cascade. 

VI. EFFECT OF A PROMOTING MODE 

In this section, we study the effect of a heat bath on 
tunneling transition of two interacting particles. In many 
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FIG. 11: Numerical solutions of transcendental equation 
m addition to the studied solution n = n, with re- 
spect to 0, there are the additional solutions, n 7^ T2, shown 
in panels (a), (b) and (c), that correspond to one, four, and 
six (pairs) additional solutions, respectively. 



tunneling reactions, interaction with a vibrational sub- 
system can often be approximated by interaction with 
a single vibrational mode (a so-called promoting mode). 
As it follows from Eq. IJ12JI . a heat bath affects only the 
dynamics of the center of mass (q± = 52)- Therefore, 
in the case of an antiparallcl motion, the medium does 
not affect the rate constant, while for parallel tunneling it 
essentially contributes to the transfer rate. For both par- 
allel and antiparallcl transfer, bilinear interaction of the 
particles with a single oscillator can make a qualitative 



change in the character of tunneling. 

At small values of the interaction parameter between 
the two tunneling particles dissipative effects become 
important. 27 In two dimensions the dissipative effects 
are more pronounced for parallel rather than antipar- 
allel tunneling. The latter increases with temperature 
with a considerable contribution to the prefactor. In the 
present work, we are interested in the tunneling rate as- 
suming only exponential evolution of the transition prob- 
ability. However, nonexponential evolution can occur in 
an nonequilibrium environment £&2£i2L2!L Such a case is 
not discussed here. Accordingly, a reservoir is assumed 
to be in thermodynamic equilibrium, i.e., the tunneling 
transition is rather slow compared to the thermodynamic 
relaxation. Thus, we assume that dissipation affects the 
value of instanton action only. 

For the case of the antiparallel tunnel transfer, the 
action H25|) can be calculated with the following vibronic 
Green's function: 



D(v n ) 



C 2 



(41) 



where lol is the frequency of the vibrational mode. Af- 
ter some tedious calculations, one obtains the following 
expression for the instanton action (|25ll : 



S = 



2u; 4 (a 2 - 6 2 )r 2uj i (a + b) 2 



O 2 
"0 







E 



/3sinh(f! (/?/2 - T- ))sinh(r> T- ) 
2fi3 s inh(O /3/2) 

0(Q 2 - Lu 2 L )cosh(n t (f3/2 - r ))cosh(fi i r ) 
(-l) 4 20f(ft 2 - ft 2 )sinh(a/?/2) 

fie fie /3sinh(ef2o) 



4(w 2 -2C 2 /lu 2 ) 
f3(n 2 - uj 2 )smh(en t 



in 2 



/3cosh(ef2o) 



^ (-l)Hftf(ft 2 - il 2 ) 4ft3sinh(/3rV2) 
x [cosh(f2 (/3/2 ~ 2t )) - cosh(O /3/ 2 )] 

A P(n 2 - ^)cosh(enQ 
~ hi ( _1 ) i40 i(°i _ ^i)smh(/3Cli/2) 
x [cosh(fi i ( / 9/2 - 2t )) + cosh(fi^/2)] } . (42) 
Here, we have introduced the following notation: 

, nl = l^+col+^^-uir + sc 2 ), 
h^+ul-y/^-uir + sc* ). 



ft; 



2 2 

bj — a 



For the parallel tunneling transition, a corresponding 
action can be found in a similar way. The action Q2J 
as a function of the parameters e* = t\ — T2 and r* = 
(ti + T2)/2, with the vibronic frequency lul and coupling 
constant C, yields (see also Ref . l27j) : 



S=(a + b){3a-b)u Z T 



u A {a + b) 2 e* Auj 2 {a+bf 
2{uj 2 -2a) fi 
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4w 4 (o + bfe* 2 to 2 {a + bf A (-l) l (w 2 - x 3 -i) 



(lu 2 - 2a)P 



2n 



coth(^) - sinh- 1 !^) (cosh[(§ - 2r*)^ 
-cosh[(^ - e*)V^i] + ^ cos h[(^ - 2t* - e*)V^i] 



i C osh[(|-2T*+ £ *V^ 



'2 



c/(a + fr) 2 
+ 2(w 2 -2a) 3 / 2 

v 2 



v 2 



'2 





-2a) 




-2a] 



+ ^cosh[(^ - 2r* - e*)V^ 2 -2a] 



1 /? 
-cosh (- - 2t* 

2 LV 2 



(43) 



Here, we have denoted 
1 



C 2 



xi, 2 = ~(w +uj l + —)t ^7, 
z u;^ z 

/ C* 2 
= J(w 2 +w£ + -2-) 2 -4w 2 t4. 



For particular values of the interaction constant a and 
in the absence of interaction with the oscillator bath, the 
critical temperature T c (at which the synchronous and 
asynchronous tunnel regimes interchange) is found from 
Eqs. i|31|l and i|38|) . These equations can be generalized 
to non-zero interaction with the promoting mode. Typi- 
cally the critical temperature is found to be in the range 
from 10 K to 400 K. In glasses, T c can be very small 
while for chemical reactions it can be rather large. Ad- 
ditionally, T c depends on the mean distance between the 
particles and, therefore, on their concentration. 

Quantum tunneling is important — when ksT c / (Tilu) < 
1. Therefore, the symmetry breaking effects can take 
place at relatively high temperatures depending on the 
"frequency" of a barrier. For example, for porphyrins the 
critical temperature T c is estimated to be about 200 K. 



ized by the different adiabatic potential energy surfaces, 
©, an( i ©i an( i made a detailed comparative anal- 
ysis of the tunneling rate for two interacting particles 
moving in parallel or antiparallel within a dissipative en- 
vironment. We have also assumed exponential evolution 
of the transition probabilityj 27 i 35 i 36 i 37 i 38 

We have shown that the change in a tunneling regime 
from synchronous to asynchronous transfer for parallel 
transition occurs as a step process, similar to phase tran- 
sition of a first order, while for the antiparallel transfer 
it resembles second order phase transition. 

We have explained the effect of a cleavage in the ex- 
perimentally observedSii 2 ^* 2 ^ 3 . temperature dependence 
of the reaction rate for two tunneling particles. It has 
been shown that the effect of symmetry breaking is sta- 
ble for the parallel and unstable for the antiparallel 
transfer, as it is observed experimentally for some por- 
phyrin systems i 2 ^ 2 ^* 2 ^ 3 . We have found a complicated 
fine structure in the bifurcation region due to the quan- 
tum fluctuations for the parallel two-dimensional tunnel- 
ing. For the antiparallel tunneling, the contribution of 
four, six, twelve, etc., pairs of trajectory becomes impor- 
tant resembling the transition to chaos. 

Additionally, we have studied interaction of two parti- 
cles with phonons. Such coupling essentially modifies the 
antiparallel and parallel transitions in different ways. As 
it follows from Eq. (|12|l . the interaction with the reser- 
voir does not change the dynamics of the center of mass 
for the antiparallel motion, while it makes a significant 
contribution to the transfer rate for the parallel transfer. 
Finally, Eq. (|31|) determines the validity condition for 
temperatures beyond of which stable two-dimensional 
synchronous tunneling correlations of all the kinds occur. 



VII. 



CONCLUSIONS 
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